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The NA62 experiment at CERN collected a large sample of K + — ¥ e + v decays during a 
dedicated run in 2007, aiming at a precise test of lepton universality by measurement of the 
^ ■ helicity suppressed ratio Rk = T(K + — > e + v)/T(K + — > ^ + v). The preliminary result of the 

D ' analysis of a partial data sample of 51089 K + — > e + v candidates is Rk = (2.500±0.016) x 10~ 5 , 

(~s , which is consistent with the Standard Model expectation. 

Introduction 

^>- ■ Due to the V — A structure of the weak interactions, the Standard Model (SM) rates of the 

leptonic meson decays P + — > £ + u are helicity suppressed. Within the two-Higgs doublet models 
(2HDM), which is a wide class of models including the minimum supersymmetric (SUSY) one, 
the charged Higgs boson (H^) exchange induces a tree- level contribution to (semi)leptonic 
decays proportional to the Yukawa couplings of quarks and leptons^M In P + — > £ + u decays, the 
H exchange can compete with the W exchange, thanks to the above suppression. 

At tree level, the H^ 1 exchange contribution to P + — > £ + u decay widths (with P = n,K, B) 
is lepton flavour independent, and is approximately described by^ 

>< : AT(P+^e+u) ^ n JMp\ 2 tan 2 /3 

%'. T SM (P+^l+v)~ \M H ) l + e tan/3' U 

Here Mjj is the charged Higgs boson mass, tan/3 is the ratio of vacuum expectation values 
of the two Higgs doublets, a fundamental parameter controlling the charged Higgs couplings, 
and £q ~ 10 -2 is an effective coupling. For a reasonable choice of the parameters (tan/3 = 40, 
Mjj = 500 GeV/c 2 ), one expects ~ 30% relative suppression of B + — > £ + v decays, and ~ 0.3% 
suppression of K + — > £ + v decays. However the searches for new physics in these decay rates are 
hindered by the uncertainties of their SM predictions. 

On the other hand, the ratio of kaon leptonic decay rates Rk = T(K e 2)/T(K^2), where 
the notation K12 is adopted for the K + — > £ + v decays, has been calculated with an excellent 



accuracy within the SM^: 
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Rf 1 = {m e /mrf "*f -« I (1 + SRqed) = (2.477 ± 0.001) x 10^. (2) 

\ m K m ii ) 

Here SRqed = (— 3.78±0.04)% is a correction due to the inner bremsstrahlung (IB) Ko,^ process. 
The ratio Rk is sensitive to lepton flavour u nive rsality violation (LFV) effe cts originating at 
one- loop level from H exchange in 2HDM , and the mixing effects in the right-handed 



slepton sector, providing a unique probe into this aspect of super symmetric flavour physics^. 
Rk receives the following leading-order contribution due to LFV coupling of the Higgs boson: 

4fe=f^Vf^ 2 |Ag|WA (3) 



i?| M \M H J \m, 

where |A^-| is the mixing parameter between the superpartners of the right-handed leptons, 
which can reach ~ 1CP 3 . This can enhance Rk by 0(1%) relative, with no contradiction to 
presently known experimental constraints (including upper bounds on the LFV r —> eX decays 
with X = 77, 7, /z//). 

The current world average (including only final results, and thus ignoring the preliminary 
NA48/2 ones) is R]^ A = (2.490 ± 0.030) x 10~ 5 , dominated by a recent measurement by the 
KLOE collaboration ^ The NA62 experiment at CERN collected a dedicated data sample in 
2007-08, aiming at a measurement of Rk with a 0.4% precision. The preliminary result obtained 
with a partial data sample is presented here. 

1 Beam, detector and data taking 

The beam line and setup of the NA48/2 experiment^ were used for the NA62 2007-08 data 
taking. Experimental conditions and trigger logic were optimized for the K^jK^a measurement. 

The beam line is capable of delivering simultaneous unseparated K + and K~ beams derived 
from 400 GeV/c primary protons extracted from the CERN SPS. Most of the data, including 
the sample used for the present analysis, were collected with the K + beam only, as the muon 
sweeping system provides better suppression of the positive beam halo component. A narrow 
momentum band of (74.0 ± 1.6) GeV/c was used to minimize the corresponding contribution to 
resolution in kinematical variables. 

The fiducial decay region is contained in a 114 m long cylindrical vacuum tank. With 
1.8 x 10 12 primary protons incident on the target per SPS pulse of 4.8 s duration, the beam 
flux at the entrance to the decay volume is 2.5 x 10 7 particles per pulse. The fractions of K + , 
7r + , p, e + and fi + in the beam are 0.05, 0.63, 0.21, 0.10 and 0.01, respectively. The fraction of 
beam kaons decaying in the vacuum tank at nominal momentum is 18%. The transverse size 
of the beam within the decay volume is 5x = 5y = 7 mm (rms), and its angular divergence is 
negligible. 

Among the subdetectors located downstream the decay volume, a magnetic spectrometer, 
a plastic scintillator hodoscope (HOD) and a liquid krypton electromagnetic calorimeter (LKr) 
are principal for the measurement. The spectrometer, used to detect charged products of kaon 
decays, is composed of four drift chambers (DCHs) and a dipole magnet. The HOD producing 
fast trigger signals consists of two planes of strip-shaped counters. The LKr, used for particle 
identification and as a veto, is an almost homogeneous ionization chamber, 27Ao deep, segmented 
transversally into 13,248 cells (2x2 cm 2 each), and with no longitudinal segmentation. A beam 
pipe traversing the centres of the detectors allows undecayed beam particles and muons from 
decays of beam pions to continue their path in vacuum. 

A minimum bias trigger configuration is employed, resulting in high efficiency with relatively 
low purity. The K e 2 trigger condition consists of coincidence of hits in the HOD planes (the so 
called Q\ signal) with 10 GeV LKr energy deposition. The K^ trigger condition consists of the 
Qi signal alone downscaled by a factor of 150. Loose lower and upper limits on DCH activity 
are also applied. 

The main data taking took place during four months starting in June 2007. Two additional 
weeks of data taking allocated in September 2008 were used to collect special data samples for 
studies of systematic effects. The present analysis is based on ~ 40% of the data sample. 



2 Analysis strategy and event selection 

The analysis strategy is based on counting the numbers of reconstructed K e2 and K^ candidates 
collected concurrently. Consequently the result does not rely on kaon flux measurement, and 
several systematic effects (e.g. due to reconstruction and trigger efficiencies, time-dependent 
effects) cancel to first order. 

To take into account the significant dependence of signal acceptance and background level 
on lepton momentum, the measurement is performed independently in bins of this observable: 
10 bins covering a lepton momentum range of [15; 65] GeV/c are used. The ratio Rk in each 
bin is computed as 

1_ N(K e2 )-N B (K e2 ) A(K^) /„ x e(gga) J_ U) 

K D N(K^ 2 )-N B (K^ 2 ) A(K e2 ) f e *e(K e2 ) / LKr ' U 

where N(K£ 2 ) are the numbers of selected Kg 2 candidates (£ = e, /i), N-q{K^) are numbers of 
background events, A(K^ 2 )/A(K e2 ) is the geometric acceptance correction, fi are efficiencies of 
e//x identification, e(Ki 2 ) are trigger efficiencies, /LKr is the global efficiency of the LKr readout, 
and D = 150 is the downscaling factor of the K^ 2 trigger. 

A detailed Monte Carlo (MC) simulation including beam line optics, full detector geometry 
and material description, stray magnetic fields, local inefficiencies of DCH wires, and time vari- 
ations of the above throughout the running period, is used to evaluate the acceptance correction 
A(K „2) / ' A(K e2 ) and the geometric parts of the acceptances for background processes entering 
the computation of Nb{K%2). The ^2(7) processes are simulated in one-photon approximation^* 
the resummation of leading logarithms^' is neglected at this stage. Simulations are used to a 
limited extent only: particle identification, trigger and readout efficiencies are measured directly. 

Due to topological similarity of K e2 and K^ decays, a large part of the selection conditions 
is common for both decays: (1) exactly one reconstructed particle of positive electric charge; 
(2) its momentum 15 GeV/c < p < 65 GeV/c (the lower limit is due to the 10 GeV LKr energy 
deposit trigger requirement in K e2 trigger); (3) extrapolated track impact points in DCH, LKr 
and HOD are within their geometrical acceptances; (4) no LKr energy deposition clusters with 
energy E > 2 GeV and not associated to the track to suppress background from other kaon 
decays; (5) distance between the charged track and the nominal kaon beam axis CDA < 3 cm, 
decay vertex longitudinal position within the nominal decay volume (the latter condition is 
optimized in each lepton momentum bin). 

The following two principal selection criteria are different for the K e2 and K^ 2 decays. K^ 2 
kinematic identification is based on the reconstructed squared missing mass assuming the track 
to be a positron or a muon: M^ im (£) = (Pk — Pe) 2 , where Pk and Pi (£ = e, //) are the four- 
momenta of the kaon (average beam momentum assumed) and the lepton (positron or muon 
mass assumed). A selection condition |M^ iss (e)| < Mq is applied to select K e2 candidates, and 
|M^ iss (/i)| < Mq for K^ 2 ones, where Mq varies from 0.009 to 0.013 (GeV/c 2 ) 2 among lepton 
momentum bins depending on M^ iss resolution. Particle identification is based on the ratio E/p 
of track energy deposit in the LKr calorimeter to its momentum measured by the spectrometer. 
Particles with 0.95 < E/p < 1.1 (E/p < 0.85) are identified as positrons (muons). 

3 Backgrounds 

K^2 decay with a mis-identified muon is the main background source in the K e2 sample. 
Sufficient kinematic separation of K e2 and K^ 2 decays is not achievable at high lepton momentum 
(p > 30 GeV/c), as shown in Fig. [T^,. The probability of muon identification as positron in that 
momentum range (E/p > 0.95 due to 'catastrophic' bremsstrahlung in or in front of the LKr) 
is P(fi -> e) ~ 4 x 10~ 6 , which is non- negligible compared to Rj^ 1 = 2A77 x 10~ 5 . A direct 
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Figure 1: (a) Missing mass squared in positron hypothesis M miss (e) vs lepton momentum for reconstructed K e 2 

and K^2 decays: kinematic separation of K e 2 and K^ decays is possible at low lepton momentum only, (b) 

Measured and simulated probability of muon identification as electron/positron P(/i — » e) vs its momentum: data 

with the Pb wall, MC simulations with and without the Pb wall (the signal region is marked with arrows). 



measurement of P(/i — > e) to ~ 10 -2 relative precision is necessary to validate the theoretical 
calculation of the bremsstrahlung cross-section^^ in the high 7 energy range used to evaluate 
the K^2 background. 

The available muon samples are typically affected by ~ 10~ 4 electron/positron contamination 
due to fi — > e decays in flight, which obstructs the P{[i — > e) measurements. In order to obtain 
sufficiently pure muon samples, a 9.2A"o thick lead (Pb) wall covering ~ 20% of the geometric 
acceptance was installed in front of the LKr calorimeter (between the two HOD planes) during 
a period of the data taking. In the samples of tracks traversing the Pb and having E/p > 0.95, 
the electron component is suppressed to a level of ~ 10 -7 by energy losses in Pb. 

The momentum dependence of P(/x — > e) for muons traversing the Pb has been measured 
with a data sample collected during a special muon (// ) run of 20h duration, and compared to 
the results of a dedicated Geant4-based MC simulation of the region downstream the spectrom- 
eter involving standard energy loss processes and bremsstrahlung^^. The data/MC comparison 
(Fig. [lb) shows good agreement in a wide momentum range within statistical errors, which val- 
idates the cross-section calculation at the corresponding precision level. The simulation shows 
that the Pb wall modifies P{[i — > e) via two principal mechanisms: 1) muon energy loss in the 
Pb by ionization decreasing P(/j, — > e) and dominating at low momentum; 2) bremsstrahlung in 
Pb increasing P(/i — > e) and dominating at high momentum. 

To estimate the K^ background contamination, the kinematic suppression factor is com- 
puted with the standard setup simulation, while the validated simulation of muon interaction 
in the LKr is employed to account for P(/i — >• e) suppression. Uncertainty of the background 
estimate is due to the limited size of the data sample used to validate the simulation. 

K(i2 decay followed by /x — >■ e decay contributes significantly to the background. How- 
ever energetic forward daughter positrons compatible to K e 2 topology are suppressed due to 
muon polarization^^. 

Ke.2-1 (SD) decay, a background by Rk definition, has a rate similar to that of K e 2~. the 
world average is BR = (1.52 ± 0.23) x 10 -5 . Theoretical rate calculations depend on the form 
factor model, and have a similar precision. Energetic positrons (E* > 230 MeV in K + frame) 
with 7 escaping detector acceptance contribute to the background. MC background estimation 
has a 15% uncertainty due to limited knowledge of the process. A recent measurement by 
KLOE , published after announcement of the NA62 preliminary result, is not used. 

Beam halo background in the K e 2 sample induced by halo muons (undergoing ^ — > e 
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Figure 2: (a) Reconstructed squared missing mass distribution M^ iss (e) for the K e 2 candidates: data (dots) 
presented as sum of MC signal and background contributions (filled areas), (b) Numbers of K B 2 candidates and 

background events in lepton momentum bins. 



Source N B /N tot 


Source N B /N tot 


Source 


N B /N tot 


K^ (6.28 ±0.17)% 
Kf,2 (A* -+ e) (0.23 ± 0.01)% 


K e2l (SD) (1.02 ±0.15)% 
Beam halo (1.45 ± 0.04)% 


K e3 
K 2 -n 


0.03% 
0.03% 


Total background: (8.03 ± 0.23)% 



Table 1: Summary of the background sources in the K e 2 sample. 

decay in flight or mis-identified) is measured directly by reconstructing K^ 2 candidates from a 
control K~ data sample collected with the K + beam dumped. Background rate and kinematical 
distribution are qualitatively reproduced by a halo simulation. The uncertainty is due to the 
limited size of the K~ sample. Beam halo is the only significant background source in the K^ 
sample, measured to be 0.25% (with a negligible uncertainty) with the same technique as for 
K e 2 decays. 

The number of K(2 candidates is N(K e2 ) = 51,089 (about four times the statistics col- 
lected by KLOEEI) and N(K^ 2 ) = 15.56 x 10 6 . The M^ iss (e) distributions of data events 
and backgrounds are presented in Fig. [2k. Backgrounds integrated over lepton momentum are 
summarized in Tabled! their distributions over lepton momentum are presented in Fig. [2Jd. 



4 Systematic uncertainties 

Positron identification efficiency is measured directly as a function of momentum and LKr 
impact point using pure samples of electrons and positrons obtained by kinematic selection of 
K + — > TT°e + v decays collected concurrently with the K e2 sample, and K\ — > ir^e^v decays 
from a special K® run of 15 hours duration. The K + and K® measurements are in good agree- 
ment. The measured f e averaged over the K e2 sample is (99.20 ± 0.05)%. Muon identification 
inefficiency is negligible. 

The geometric acceptance correction AiK^)! A{K e2 ) is strongly affected by the radia- 
tive K e2l (IB) decays. A conservative systematic uncertainty is attributed to approximations 
used in the K e2 y IB simulation. The resummation of leading logarithms'^ is not taken into ac- 
count, however no systematic error is ascribed due to that. An additional systematic uncertainty 
reflects the precision of beam line and apparatus description in the MC simulation. 

Trigger efficiency correction e{K e2 )/e{K^ 2 ) = 99.9% accounts for the fact that K e2 and 
K(j,2 decay modes are collected with different trigger conditions: the E > 10 GeV LKr energy 
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Figure 3: (a) Measurements of Kk in lepton momentum bins, (b) The world average of Rk- 
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Statistical 
K^ 2 



0.012 
0.004 



K e2l (SD) 0.004 



Beam halo 
Positron ID 
IB simulation 



0.001 
0.001 
0.007 



Geom. acceptance 0.002 

Trigger dead time 0.007 



Table 2: Summary of uncertainties of Rk'- statistical and systematic contributions. 

deposition signal enters the K e2 trigger only. A conservative systematic uncertainty of 0.3% is 
ascribed due to effects of trigger dead time which affect the two modes differently. LKr global 
readout efficiency /LKr is measured directly to be (99.80 ± 0.01)% and stable in time using 
an independent LKr readout system. 

5 Result and conclusions 

The independent measurements of Rk in lepton momentum bins, and the result combined 
over the momentum bins are presented in Fig. [3h- The uncertainties of the combined Rk are 
summarised in Tabled The preliminary result is Rk = (2.500 ±0.012 sta t. ±0.011 syst .) x 10~ 5 = 
(2.500 ± 0.016) x 10~ 5 , which is consistent with the SM expectation. Analysis of the whole 
2007-08 data sample is expected to decrease the uncertainty of Rk down to 0.4%. A summary 
of Rk measurements is presented in Fig. [3b: the current world average is (2.498 ±0.014) x 10 -5 . 
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